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ABSTRACT 

We investigate the magnetic field evolution in a series of galaxy minor mergers using 
the N-body /smoothed particle hydrodynamics (SPH) code Gadget. The simulations 
include the effects of radiative cooling, star formation and supernova feedback. Mag- 
netohydrodynamics (MHD) is implemented using the SPH method. We present 32 
simulations of binary mergers of disc galaxies with mass ratios of 2:1 up to 100:1, 
whereby we have additionally varied the initial magnetic field strengths, disc orienta- 
tions and resolutions. We investigate the amplification of a given initial magnetic field 
within the galaxies and an ambient intergalactic medium (IGM) during the interac- 
tion. We find that the magnetic field strengths of merger remnants with mass ratios 
up to 10:1 saturate at a common value of several /iG. For higher mass ratios, the field 
strength saturates at lower values. The saturation values correspond to the equipar- 
tition value of magnetic and turbulent energy density. The initial magnetization, disc 
orientation and numerical resolution show only minor effects on the saturation value 
of the magnetic field. We demonstrate that a higher impact energy of the progeni- 
tor galaxies leads to a more efficient magnetic field amplification. The magnetic and 
turbulent energy densities are higher for larger companion galaxies, consistent with 
the higher impact energy supplied to the system. We present a detailed study of the 
evolution of the temperature and the bolometric X-ray luminosity within the merg- 
ing systems. Thereby we find that magnetic fields cause a more efficient increase of 
the IGM temperature and the corresponding IGM X-ray luminosity after the first 
encounter. However, the presence of magnetic fields does not enhance the total X-ray 
luminosity. Generally, the final value of the X-ray luminosity is even clearly lower for 
higher initial magnetic fields. 

Key words: methods: N-body simulations - galaxies: spiral - galaxies: magnetic fields 
- galaxies: evolution - galaxies: kinematics and dynamics 



1 INTRODUCTION 

In the framework of hierarchical growth of structure in the 
universe, galaxy interactions are believed to be an essen- 
tial part of galaxy formation and evolution. The Lambda 
cold dark matter (ACDM) cosmology predicts the forma- 
tion of dark matter haloes due to gravitational instabili- 
ties in the early universe which later on form larger haloes 
by gas accretion an d halo mergers jWhite fc Reesl 1 19781 : 
IWhite fe Frenklll99ll ). In the continuing process, baryonic 
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particles get gravitationally bound to the haloes, forming 
structures which are consistent with the galaxies observed 
today. The merger rate of dark matter haloes is an increas- 
ing fu nction of redshift jKolatt et al.|[l999l : iGottlober et al.l 
l200l[ ). implying that collisions and mergers were much more 
frequent in the earl y Universe. Moreov er, simulations of 
merger history trees ( Wechsler et al.|[2002t ) indicate that ma- 
jor merger events were comparatively rare, and thus dark 
matter haloes grow mainly by the accretion of smaller ob- 
jects. More precisely, minor mergers are expected to be at 
least on e order of magnitude mo re common than major 
mergers l|Hernauist fc Mihoslll995l . and references therein). 
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The most obvious consequence of galaxy collisions, the 
shaping of the part i cipati ng objects, was first described by 
iTbomre fc Toomra (|l972h . Interactions of galaxies lead to 
significant changes of the dynamics of the progenitor galax- 
ies due to the alteration of the gravitational pot ential (e.g. 
iToomre fc Toomrd Il97i iNaab fc Burkertl l2003h . Thereby, 
most of the galaxy collisions result in a merger of the pro- 
genitor systems. 

So far, simulations of galaxy mergers were predominantly 
dedicated to studies of star formation, stellar dynamics, 
gas flows, supermassive black hole s or fe e dback from stars 
and bl ack holes (e.g . ICox et al] 120081: [pTlVlatteo et al l 



20051; 



.Springe! et al l2005al lbl; | on et alj 

Johansson ct al ] 12003 ). However, galaxy mergers are also 
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interesting in the context of the amplification and restruc- 
turing of small-scale magnetic fields within t he scope of the 
global evol ution of cosmic magnetism (e.g. iKotarba et al.l 
I2009l . l2010h . Observations of galaxies have revealed that vir- 
tually all galaxies host magne tic fields, with fi eld strengths 
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between I and 10 /xG (see e.j 
IChvzv et al1l2007l : IVollmer et al]|2010l ). 

For comparison, the intergalactic ma gnetic field is usually es - 
timated to be less than 10~ 8 G (e.g. iKronberg et al.|[20ol l. 
The comparatively strong galactic magnetic fields are com- 
monly expla ined by the act ion of a galactic dynamo (for a re- 
view see e.g. lKulsrudlll999l and references therein). This dy- 
namo, which is derived using the mean field theory, acts on 
the basis of the conversion of poloidal into toroidal fields and 
vice versa. Thereby, turbulent motions generated by stellar 
activity are transformed into cyclonic motions due to Corio- 
lis forces. The cyclonic motions are responsible for the con- 
version of toroidal flux into poloidal flux (a-effect). The dif- 
ferential rotation of the galactic disc causes the conversion of 
poloidal into toroidal flux (fi-effect). This afi dynamo leads 
to an efficient amplification of the magnetic field. A simi - 
lar approach exploits the Parker instability l|ParkeilH992l ') 
caused by cosmic ray (CR) particles to generate large-scale 
turbulence which in turn results in an efficient dynamo ac- 
tion. This CR-driven dynamo process has been recently sim- 
ulated bv lHanasz et al.l |2009l ). showing an equipartition of 
the turbulent energy with magnetic field energy and cosmic 
ray energy within the gas. However, despite their success 
in explaining the magnetic field strengths and structures in 
fully evolved local galaxies, these dynamos are challenged by 
recent ob servations of strong (fiG) magnetic fields at high 
redshifts (|Bernet et al.|[20o3 ). Therefore, further amplifica- 
tion processes of galactic magneti c fields acting on short 
timescales have to be investigated. lArshakian et al.l (|2009T ) 
used analytical considerations to show the importance of 
turbulence in protogalactic haloes, which may lead to ef- 
ficient small-scale dynamo action. This work thus already 
showed the importance of turbulence for the evolution of 
magnetic fields in the early universe. H owever, the conclu- 
sions made by lArshakian et al.1 (|2009h are based on only 
rough estimates of the amplification timescales and the tur- 
bulence driven by galactic collisions. 

In order to investigate the idea of an i nteraction-driv e n am - 
plification of galactic magnetic fields, IKotarba et"ai1 (|2010h 
have performed numerical simulation of a galactic (major) 
merger of the Antennae Galaxies, including the evolution of 
magnetic fields. They found that the magnetic field within 
the colliding system gets amplified by compression and shear 



flows up to a saturation value of ~ 10/iG, independent of 
the initial magnetic field of the progenitor discs, which was 
varied between 10 -9 and 10~ 6 G. Within this work, the sat- 
uration level was found to be near equipartition between 
magnetic and turbulent gas pressure, in good agreement 
with theoretical expec tations of the turbulent dynamo (e.g. 
lArshakian et al. 20091. and refer ences therein). In a continu- 
ative study, Kotarba et al.|[2~01ll considered a major collision 
of three galaxies. These studies confirmed the saturation of 
the galactic magnetic field at the equipartition level of sev- 
eral fiG independent of the initial magnetic field. Further- 
more, an additionally included ambient IGM allowed also 
for studies of its magnetization, which was shown to sat- 
urate at « 10 _8 G. However, all of these studies are dedi- 
cated to galactic major mergers. As galaxy minor mergers 
are expected to be far more frequent within the process of 
structure formation, it is definitely interesting to consider 
the influence of the mass ratios of the progenitor galaxies on 
the magnetic field amplification and saturation value. This 
idea is pursued within the presented work. 
Moreover, IKotarba et~ai]|201ll also found that the initial 
magnetization of the galaxies and the IGM affects the 
propagation of interaction-driven shocks within the IGM. 
Thereby, the shocks are stronger (gaining higher Mach num- 
bers) for higher initial magnetizations. This is also indicated 
by a larger extent of the shock heated regions and the higher 
temperatures of these regions. This hot gas is expected to 
radiate in X-rays due to thermal bremsstrahlung. 
Observations w i th the ROSAT X-ray telescope 
|Read fc Ponmanl 19981) and the Chandra X-ray obser- 
vatory (e.g. iFabbiano et al.l l200lf ) of interacting galaxy 
systems revealed the production of extended X-ray emission 
during the interactions. Simulations focusing on the X-ray 
emission o f inter acting galaxies were recently performed b y 
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2006) 1 and ISinha fc Hollev-Bockelmannl l|2009l ). 
( 20061 ') considered the X-ray emission generated 
by interactions of galaxies including gaseous discs, but 
without including an ambient IGM gas. They found that 
the cold gas in the galactic discs gets shock-heated to 
X-ray emitting temperatures during the collision. The 
X-ray emission in their simulations (different major merger 
scenarios with different gas fractions) increases during 
the progress of the collision and peaks at the time of the 
final merger. They showed that galaxy mergers are able 
to generate X-ray haloes of elliptical merger remnants. 
In contrast to this investigations of cold gaseous discs, 
ISinha fc Hollev-Bockelmannl l|2009f ) simulated the collisions 
of galaxies including hot gas within the galactic haloes. 
However, their models did not contain any disc gas. In their 
simulations of different major and minor merger (3:1 and 
10:1 mass ratio) scenarios with different gas fractions and 
orbital properties, the gas in the haloes gets shock-heated 
already before the first encounter. The strongest shocks and 
highest X-ray emissions occur after the first encounter of 
the galaxies. They also find that the peak X-ray luminosity 
increases with increasing progenitor masses and also with 
increasing gas fraction. 

Since the increase of the IGM temp erature also depends on 
the initial magnetic field strength |Kotarba et al.l [201 lh ■ a 
connection between the X-ray emission and the magnetiza- 
tion may be expected, but detailed investigations concerning 
this interconnection are still missing. Also, there are no stud- 
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ies of the X-ray emission during galactic interactions which 
include both gaseous discs and a sourrounding hot IGM gas. 
Therefore, within the presented work we want to address 
the following issues: Until now, simulations of galactic merg- 
ers including magnetic fields have only been performed for 
major merger scenarios. However, as minor mergers are ex- 
pected to be the more frequent events, we aim to study 
the effect of the mass ratio of the progenitor galaxies on 
the magnetic field evolution. Based on 32 different mi- 
nor merger scenarios, we investigate the amplification of 
a given initial magnetic field within the galaxies and an 
ambi ent IGM numerically usi ng the MHD version of Gad- 
get (|Dolag fc Stasvszvnl200d ). A particular focus is thereby 
placed on the evolution of the expected X-ray emission of the 
simulated systems. We also discuss the resolution-dependent 
numerical divergence in our simulations. 
The paper is organized as follows: In section 2 we briefly de- 
scribe the numerical method. The galaxy models and merger 
scenarios are described in detail in section 3. In section 4 we 
present the results of our simulated merger scenarios, par- 
ticularly the magnetic field evolution and X-ray emission. 
Finally, we summerize and discuss our findings in section 5. 



2 METHOD 

All galaxy collision simulations ha ve been performed 
with the N-bo dy/SPH code Gadget jSpringel et al.ll200ll : 
ISpringel|[2005l ').' This code models dark matter and stars 
as a self-gravitating, collisionless fluid, which is treated 
with the traditional N-body approach. Gravitational in- 
teractions are computed with a Barnes & Hut tree con- 
struction jBarnes fc HutJ 119861 ). The IGM and interstellar 
medium (ISM) gas is described as a conductive, ideal fluid, 
whereby hydrodynamics a re treated with the SPH method 
(for recent reviews see e.g. iRosswog S.ll2009l : ISpringel|[2O10l : 
IPricd l201lh . The ap plied SPH formulation conserves both 
energy and entropy ijSpringel fc Hernauistll2002h. An addi - 
tional MHD implementation by iDolag fc Stasvszvnl (|2009l ) 
allows for the evolution of magnetic fields. Below, we briefly 
describe the SPH and smoothed particle magnetohydrody- 
namics (SPMHD) method s. Fo r detai l ed des criptions see e.g. 
IDolag fc Stasvszvnl (|2009i > and lPricd ^Olll ). 



2.1 SPH 

In SPH, the fluid is discretized by mass elements rrii at po- 
sitions Xi. The density can then be estimated using 



Pi 



(1) 



where the kernel W jMonaghan fc Lattan zio 1985) is a 
weighting function determining the contribution of the par- 
ticles inside a smoothing sphere with the radius hi: 
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To enable spatial adaptivity, the smoothing length hi is eval- 
uated iteratively together with the density via 

l/d 
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whereby the parameter 77 is defined by the dimension d 
of the problem and the shape of the kernel, within our 3- 
dimensional simulations by a cubic kernel: 
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with the number of neighbouring particles N. 

The equation of motion for the SPH particles can be derived 

using a variational principle, which leads to 



dvi 



hyd 



p\ 



(5) 



with Wt — W(xi — Xj,hi) and Wj = W(xi — Xj,hj). The 
correction terms ff° account for the adaptive smoothing 
lengths: 



/CO 
i 



^ K dpi 
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We use the standard implementation of the artificial viscos- 
ity with values of a — 2 a nd ft = 1.5 for the dimensionless 
parameters (for details see |Price]|201ll ). 



2.2 SPMHD 

In ideal MHD, the evolution of a magnetic field can be fol- 
lowed using the induction equation 



^ = (B-V)v-B(V-v), 



(7) 



whereby the constraint V ■ B = has been used. The SPH 
discretization of the induction equation reads 

=fr-Yrn j [B i (v ij -ViWi)-v ij (B i -V i W i )]. (8) 
at pi t-r 1 

3 

We also apply artificial magnetic dissipation as de- 
scribed in Dolag fc Stasvszvnl (|2009h which is based on 
IPrice fc Monaghanl l|2004al ). 

The magnetic field exerts a feedback on the plasma via the 
Lorentz force. The magnetic contribution to the acceleration 
of a particle can be written in the following symmetric and 
conservative form 

■ fc V nag 18M kl 



dv' 
~~dt 



p dx l 

with the magnetic stress tensor |Phillips fc Monaghanlll985l ) 

Po 



-B*5 kl 



(10) 



Hence, the contribution to the equation of motion 
caused by magnetic field s can be discretized as follows 
(|Price fc Monaghanl |2004bT ) 



dt 



ft 



(ii) 

However, for strong magnetic forces, this momentum con- 
serving form can lead to numerical instabilities, i.e. clumping 
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of the particles (|Phillips fe Monagharj[l985h , because of the 
non- vanishing numerical divergence of B, which will not be 
stabilized in cases where the magnetic pressure exceeds the 
gas pressure. One method to cir cumvent this problem was 
suggested by lB0rve et alj |200ll ). Within this method, the 
effects of any unphysical source terms of B are subtracted 
from the momentum equation. The corresponding correction 
term can be calculated via 



dvi 



MO ^ 



fr Ei ViWi + fr Ei. v . 

Mi Hj 



(12) 



where (3 is a constant factor, with a typical value of f3 = 1. 
Additionally, a threshold for the divergence force subtraction 
was introduced in order to account for situations in which 
the acceler ation due to the div ergence force could become 
dominant (|Kotarba et al.ll201ll . Stasyszyn & Dolag 2011, in 
preparation). This threshold is set to half of the value of the 
current Lorentz force. Generally, it was found that the ef- 
fects due to the violation of momentum conservation by the 
divergence force subtraction are of an insignificant order, 
whereas the implementation of the cor rection term signifi- 
cantl y improves the results in test cases (|Dolag fe Stasvszvnl 
l2009h . 

The MHD version of Gadget was already successfully 
employed for the analysis of magnetic fields in star for- 
mation, in isolated galaxies and me rgers of galaxies as 
well as within gala x y clusters (e.g. iBiirzle et al.l 2011bl; 
Kotarba et all 120091 ; foonnert et alTl2009l ; iBonafede et"ail 



201ll ). 



2.3 



Radiative cooling, star formation and 
supernova feedback 



Radiative cooling is modeled assuming the gas to be 
optically thin an d in collisional ionization equilibrium 
l|Katz et all 119961 ). The gas is assumed to be of primor- 
dial composition (hydrogen fraction X=0.76, helium frac- 
tion Y=0.24), hence, metal-dependent effects are neglected. 
This approach is reasonable as we will consider the primary 
mechanism of X-ray emission (see chapter 4.3) to be ther- 
mal bremsstrahlung, which is c onsistent with th e applica- 
tion of zero-metallicity cooling |Cox et al.l l2006). Also, we 
account for a homogeneous extra galactic ultraviolet (UV ) 

Hvladaulll996h . 
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background in our simulations (e.g. 
Star formation (SF) and associated supernova feedback are 
treated by the hybrid mul tiphase model as described in 
ISpringel fe HernquistJ (|2003T ). where the interstellar medium 
is modeled as a two-phase fluid (|McKee fc Ostrikerlll977t ). 
This fluid is assumed to consist of cold clouds which are 
embedded in an ambient hot medium. At high densities, 
material gets bound in these clouds and star formation can 
take place, which in turn lowers the density of the ambient 
gas phase. This leads to a decrease of its radiative losses 
and enables the heating of the region caused by supernova 
explosions, resulting in pressure support for the ISM. The 
multiphase model is realized numerically by describing each 
star-forming particle to be made up of two components, 
a hot diffuse component and a cold component, whereby 



the local gas density and temperature determine the rela- 
tive amounts of the phases. The effective equation of state 
EOS) for the gas is given by P eff = (7 - ■ l)(p h i ih + p c u c ) 
Springel fc Hernquistll2003l ; ISpringel et al.ll2005br i, with the 
density and the specific thermal energy of the hot (ph and 
ith) and the cold phase (p c and u c ), and an adiabatic index 
of 7 = 5/3. 



3 INITIAL CONDITIONS 
3.1 Galaxy models 

In order to perform a series of simulations of unequal mass 
mergers, we set up structurally similar galaxy models with 
different total masse s. The galaxies are se t up using the 
method described by ISpringel et all (|2005bT ). This method 
allows for a galaxy model consisting of a cold dark matter 
halo, an exponential stellar disc, a stellar bulge (all of these 
components being collisionless N-body particles) and an ex- 
ponential gaseous disc (SPH particles). Both the virial mass 
M200 and the virial radius r2oo of the galaxy depend on the 
virial velocity V200 (|Springel et al.l l2005bl ; Ijohansson et al.l 
I2009D via 



T200 



M200 



^200 

WHo' 



vjoo 
WGHo ' 



(13) 



(14) 



with the Hubble constant Ho=/iT00 km s _1 Mpc -1 and h — 
0.71. Thus, we determine the masses of our galaxy models 
by varying the virial velocity V200 ■ We assume the same halo 
concentration CC, spin parameter A, disc and bulge mass 
fractions md,mb, disc spin fraction gas fraction f, disc 
height zo, bulge size lb and the scale length of extended 
gas disc l s for all of the used galaxy models. These common 
parameters and the parameters of the multi-phase model are 
given in Table 1, where masses are given in units of the total 
galactic mass M to t and scale lengths in units of the stellar 
disc scale length 

We set up 12 medium resolution galaxy models, 4 high res- 
olution models and two low resolution models. Model Ml is 
the most massive galaxy in our sample, and the other mod- 
els are numbered according to their mass with respect to 
the Ml model, i.e. model M2 is a galaxy with half the mass 
of model Ml, M3 has one third of the mass of Ml and so 
forth 1 . This allows for a transparent merger denotation, i.e. 
merger M1M2 describes a 2:1 merger, M1M3 a 3:1 merger 
and so on. We refer to the mass ratios in terms of the mass 
of the companion galaxy divided by the mass of Ml. 
The used particle numbers, total masses, virial radii and 
stellar disc scale lengths used for our different galaxy mod- 
els are summarized in Table 2. These setups result in a 
mass of the gas particles of m gas ~ 4 • 10 4 h^ 1 Mq for 
the medium resolution, m gas ~ 4 ■ 10 3 Mq for the high 



1 By varying the virial velocity it is not possible to set the mass of 
the galaxies to an exact predefined value. Therefore, our smaller 
models contain only approximately a certain percentage of the 
model Ml, which is reflected in the irregular particle numbers in 
Table 2. 
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Table 2. Galaxy Model Parameters 



MEDIUM RESOLUTION 



Model 


Mtot 
[10 10 M m ] 


R.200 

[kpc/h] 


Id 

[kpc/h] 


jya 1 


^disc 


7V C 

gas 


^bul e 

u ge 


Ml 


134.14 


160.0 


7.09 


400 000 


960 000 


240 000 


400 000 


M2 


63.96 


125.0 


5.54 


200 000 


457 764 


114 441 


190 735 


M3 


43.60 


110.0 


4.87 


133 333 


311 954 


77 988 


129 981 


M4 


32.75 


100.0 


4.24 


100 000 


234 375 


58 594 


97 656 


M5 


26.77 


93.5 


4.14 


80 000 


191 579 


47 895 


79 824 


M6 


22.32 


88.0 


3.90 


66 667 


159 720 


39 930 


66 550 


M7 


19.07 


83.5 


3.70 


57 143 


136 449 


34 112 


56 854 


M8 


16.77 


80.0 


3.54 


50 000 


i on f\r\r\ 

120 000 


30 000 


50 000 


iviy 


1 A 

±4.yo 


1 t .u 


"i A^ 


AA AAA 
44 444 


1 07 nnn 
1U i uuu 


ZD 1 OU 


AA ^8*3 
44 ooo 


M10 


13.27 


74.0 


3.28 


40 000 


94 974 


23 744 


39 573 


M15 


8.99 


65.0 


2.88 


26 667 


64 365 


16 091 


26 819 


M20 


6.73 


59.0 


2.61 


20 000 


48 136 


12 034 


20 057 








HIGH 


RESOLUTION 








Ml 


134.14 


160.0 


7.09 


4 000 000 


9 600 000 


2 400 000 


4 000 000 


M3 


43.60 


110.0 


4.87 


1 333 333 


3 119 536 


779 884 


1 299 807 


M50 


1.91 


43.5 


1.93 


80 000 


192 921 


48 230 


80 384 


M100 


0.96 


34.5 


1.53 


40 000 


96 243 


24 061 


40 101 








LOW 


RESOLUTION 








Ml 


134.14 


160 


7.09 


40 000 


96 000 


24 000 


40 000 


M3 


43.60 


110 


4.87 


13 333 


31 195 


7 799 


12 998 


(a) collis 


onless particles within dark matter halo (b) collisionless particles 


within disc 




(c) gas particles within d 


sc 




(d) collisi 


onless particles within bulge 





Table 1. Parameters of initial setup common to all galaxy models 



Table 3. Gravitational softening lengths e 



Disc Parameters 



Concentration 


CC 


12 


Spin parameter 


A 


0.1 


Disk mass fraction" 


m d 


0.05 ■ Mtot 


Bulge mass fraction™ 


in]., 


0.02 • Mtot 


Disk spin fraction 


.id 


0.05 


Gas fraction 


f 


0.2 


Disk height™ 


zo 


0.2 ■ l d 


Bulge size™ 


h 


0.2 ■ l d 


Scale length of extended gas disc™ 


k 


6 ■ ; d 


Multi Phase Model Parameters 


Gas consumption timescale 


tMP 


8.4 Gyr 


Mass fraction of massive stars 


Pmp 


0.1 


Evaporation parameter 


A 


4000 


Effective supernova temperature 


TsN 


4 ■ 10 8 K 


Temperature of cold clouds 


T CC 


1000 K 



(a) Mtot and 1^ are given in Table 3. 



resolution and m gas ~ 4-10 h~ Mq for the low resolution. 
For the medium resolution, the fixed gravitational softening 
lengths are e = 0.11 /i~*kpc for the dark matter particles 
and e = 0.022 ft" 1 kpc for the gas, bulge and halo particles. 
For the low and high resolution runs, these softening lengths 
have been adjusted using h ncw — h \d • (AW/A^ew) (see 





Halo 


Disk 


Gas 


Bulge 




[pc/h] 


[pc/h] 


[pc/h] 


[pc/h] 


low resolution 


240 


48 


48 


48 


medium resolution 


110 


22 


22 


22 


high resolution 


52 


10 


10 


10 



e.g. lDehnen|[200l1 : Ijohansson et al.ll2009t ) . The gravitational 
softening lengths are summarized in Table 3. The minimum 
SPH smoothing length for the gas particles is l.Oe. 



3.2 Galactic magnetic field 

The initial magnetization of the progenitor discs is given 
by _B X = -B g ai,o an d B y — B z = G with the z-axis being 
the axis of rotation. Thus, the initial field lies always in 
the plane of the galactic discs. As we are interested in the 
influence of minor mergers on the galactic magnetic field 
evolution, particularly its amplification, we mainly focus on 
a small initial magnetic field strength of B ga i,o = 10 -9 G. 
This value is by three orders of magnitude small er than the 
typic al observed galactic magnetic field value (|Beck et al.1 
Il996t ). However, we also simulate two of our merger scenarios 
(see section 3.6.) with a initial magnetic field strength of 
B gal.O = 10~ 6 G as well as mergers without any magnetic 
fields for comparison. 
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Table 4. Simulated galaxy merger scenarios 



Merger scenarios 



Scenario 


mass ratio 


orbit 


initial -Bgai.o 


initial Bigm,o 


resolution 


M1M2 G9I9 


2:1 


prograde 


10 -9 G 


10~ 9 G 


medium 


M1M3 G9I9 


3:1 


prograde 


10 -9 G 


10~ 9 G 


medium 


TV It 1 TV It A f~i ATn 

M1M4 G919 


4:1 


prograde 


10 -9 G 


10 -9 G 


medium 


tv it 1 a it r oTn 

M1M5 G919 


5:1 


prograde 
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3.3 Orbits 



3.4 IGM 



Within all of the presented simulations the largest (i.e. most 
massive) galaxy model Ml interacts with one of the smaller 
(i.e. less massive) models M2-M100. As only minor mergers 
are considered, we will refer to the smaller galaxy in each 
simulation as the "companion galaxy". In order to ensure a 
collision, the galaxies are set on a parabolic orbit, result- 
ing in a prograde encounter (i.e. the spin direction is the 
same within both galaxies) in most of the simulations. Two 
of the models are set on a retrograde encounter for compari- 
son. The initial separation r scp of the galaxies is determined 
by the sum of their virial radii. The pericenter distance in 
all simulations is r p = 5 /i _1 kpc. The disc orientation (see 
iToomre fc Toomrej Tl973 ) is set to i = 60° and w = 60° 
for the largest galaxy Ml and i = 60° and uj — —60° for 
the companion galaxy, respectively, within all of our merger 
scenarios. A detailed description of these merger scenarios 
is given in section 3.6. 



We include an ambient IGM composed of additional gas 
partic les surrounding the galaxies similar to iKotarba et al.l 
(|201il ). The IGM particles are arranged in a hexagonal 
closed-packed lattice. The mass of the IGM gas particles is 
the same as the mass of the galactic gas particles. The vol- 
ume filled with the IGM is 700 x 700 x 700 h~' A kpc 3 centered 
on the common center of mass of the progenitor galaxies at 
the beginning of the simulation. For simplicity, we assume 
that the IGM is pervading the galaxies. The density of the 
IGM is pigm = 10 -29 g cm -3 , resulting in IGM particle 
numbers of Ni GM = 89 713 for the low, ]Vigm= 905 205 for 
the medium and Wigm = 8 997 083 for the high resolution 
setups. 

The internal energy of the IGM is calculated via uigm = 
t>2oo/2; using the virial velocity V200 °f the larger progenitor 
galaxy Ml. This sets the temperature of the IGM to the 
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Figure 1. Evolution of the projection of the mean line-of-sight density at eight different time steps (given within the panels) for the 
minor merger scenario M1M4_G6I9. 



virial temperature at the virial radius of the galaxy Ml: 



2 m pf i 

J IGM = — ^IGM— ; 

3 Kb 



1 2 N m P n 
3 Kb 



(15) 



with the mean molecular weight for a fully ionized gas of pri- 
mordial composition /i ~ 0.588, proton mass m p and Boltz- 
mann constant Kb- 



3.5 IGM magnetic field 

The initial magnetic field of the IGM is assumed to be homo- 
geneous and directed in x-direction, i.e. Bigm,o — £?igm,x, 
with now the x-y-plane being the orbital plane. Due to our 
setup, the IGM magnetic field is also pervading the galaxies. 
Bigm,o is assumed to have values of either 10~ 9 G or 1CP 12 
G, respectively, depending on the merger scenario (see the 
next section). For the runs excluding magnetic fields also 
the IGM is not magnetized. 



3.6 Merger scenarios 

We have performed 11 standard, medium-resolution minor 
merger simulations with a standard galactic and IGM mag- 
netic field of -Bgai.o — £>igm,o = 10~ 9 G and a prograde disc 
orientation. The mass ratios of the discs range from 2:1 down 
to 20:1. The denotation "G9I9" in our scenario titles refers 
to the exponent of the initial galactic (G) and IGM (I) mag- 
netic field strengths. We resimulate all of these 11 merger 
scenarios with a lower IGM magnetic field of Bigm.o = 10~ 12 
G. We also perform two additional simulations with mass ra- 
tios of 50:1 and 100:1, respectively. For these two simulation 
we have to adopt a high resolution in order to sample the 
companion galaxy with at least 10 000 gas particles. Two 
of the standard scenarios (mass ratios 4:1 and 7:1) are ad- 
ditionally resimulated with a higher galactic magnetic field 
of Bgai.o = 10~ 6 G (we will refer to these models also as 
to the present-day mergers). Also, we resimulate two of the 
standard scenarios (mass ratios 3:1 and 7:1) with a retro- 
grade disc orientation. Additionally, we resimulate the 3:1 




t [Gyr] 

Figure 2. Angular momentum of the larger progenitor galaxy 
Ml as a function of time for the standard merger scenarios with 
^gal,o = ^IGM.O = 10~ 9 G. 

standard scenario with both low and high resolution for res- 
olution studies. Finally, two of the standard scenarios (mass 
ratios 4:1 and 7:1) are resimulated excluding magnetic fields 
completely. Thus, we present altogether 32 merger scenarios 
in this paper. The different scenarios are listed in Table 4. 



4 SIMULATIONS 

We let our initial system evolve for 200 Myr to allow possible 
numerical discontinuities associated with the initial setup 
(e.g. effects caused by overlayed magnetic fields of galaxies 
and IGM) to relax. 

In the following, the physical parameters shown in the plots 
are calculated by taking the mean value (L) — Lj of 

all particles, if not stated otherwise. 



4.1 Global evolution 

Fig. 1 shows the projection of the mean line-of-sight den- 
sity exemplarily for the present-day scenario M1M4_G6I9, 
whereby the 8 panels display a sequence of increasing time. 
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In the beginning of the simulation, the galaxies are moving 
towards each other due to their mutual gravitational attrac- 
tion. The first encounter takes place at t « 1.3 Gyr 2 , where- 
upon prominent tidal arms are developing. The first en- 
counter is followed by a series of subsequent encounters until 
the final merger occurs at about t ~ 4 Gyr. For smaller pro- 
genitor galaxies, subsequent encounters take place at later 
times (and more of them until the final merger) due to the 
weaker gravitational attraction (not shown). For example, 
the final merger takes place at 3.4, 6.5, and 8.5 Gyr for the 
M1M2, the M1M7 and the M1M10 scenario, respectively. 
In most of the simulated scenarios, the disc of the larger 
progenitor galaxy Ml outlasts the interaction. However, in 
scenarios with the largest companion galaxies M2 and M3 
(scenarios M1M2 and M1M3), the disc of Ml gets largely 
disrupted. 

In order to quantify the amount of disc disruption, we esti- 
mate the angular momentum of the galaxy Ml by calculat- 
ing L = m(r x v) for all particles initially belonging to Ml, 
with the L-axis being perpendicular to the galactic plane. 
Given this simple approach we may miss particles which 
are accreted onto the galaxy during the simulations, and 
spuriously account for particles which have already left the 
disk. Hence, the calculated angular momentum cannot rep- 
resent the true angular momentum perfectly. However, its 
evolution can nevertheless reveal the dynamical evolution of 
the disc. Fig. 2 shows the normalized angular momentum 
|L|/|Lo| (with Lq the angular momentum at the beginning 
of the simulations) of Ml with time. At the time of the 
first encounter (f» 1.3 Gyr), the angular momentum gener- 
ally decreases, whereby the decrease is stronger for higher 
mass ratios. This is reasonable as larger companion galax- 
ies can disrupt the disc of galaxy Ml more heavily. For the 
mergers with the largest mass ratios, only a slight decrease 
in the angular momentum is visible. The fluctuations of the 
angular momentum within all scenarios are due to the asym- 
metric mass distribution within the gaseous discs after the 
first encounter. They are not caused by further encounters. 
These fluctuations are more pronounced for larger compan- 
ion galaxies, as the encounters have a stronger impact on 
the mass distribution. 



4.2 Magnetic field evolution 

The magnetic field is as sumed to get enhanced through tur- 
bulent motions (see e.g. iBrandenburg fc Subramanianll2005l . 
for a review) , which are caused by the impact and rotational 
energies of the galaxies. Therefore, we expect galaxy colli- 
sions with higher impact energies (i.e. galaxies with larger 
companion galaxies, see chapter 4.3.1) to be able to amplify 
the magnetic field more efficiently. Also, collisions of systems 
with different spin directions may be expected to lead to a 
more efficient amplification, as different rotational velocities 

2 All of the scenarios show the first encounter at t 1.3 Gyr, 
but the first encounter for smaller mass ratios is slightly earlier 
than for collisions with larger mass ratios. The reason for this 
is the separation distance of the galaxies, which we set to the 
sum of the virial radii of the two galaxies. As the dark matter 
haloes already overlap at this distance, the Keplerian orbit is 
only nearly-parabolic. 



can cause higher shear motions and could thus lead to a 
higher turbulence within the system. 

Fig. 3 shows the morphological evolution of the projection 
of the mean line-of-sight total magnetic field again for the 
present-day scenario M1M4_G6I9, whereby the 8 panels 
display a sequence of increasing time. Before the first en- 
counter at about 1.3 Gyr, the galactic magnetic field gets 
wound up by the differential rotation of the discs and gains 
a non-axisymmetric pattern with two magnetic arms (Fig. 
3, left upper panel). Shocks and interaction-driven outflows 
caused by the first encounter are propagating into the IGM, 
whereby the IGM magnetic field is strengthened within the 
shocked regions (Fig. 3, s econd upper panel) , a n effect which 
has been also found by iKotarba et all (|201 lh . During the 
subsequent evolution, further encounters take place, which 
are also accompanied by shocks and outflows and therefore 
lead to a further magnetization of the IGM. At the time of 
the final merger (~ 4 Gyr), the magnetic field in the galax- 
ies (i.e. within the region with a density higher than 10 -29 
g cm -3 ) has approximately retained its initial value of 10 -6 
G, whereas the IGM magnetic field got amplified within an 
extended region around the galaxies up to a value of several 
1(T 9 G. 

4-2.1 Dependence on the mass ratio 

Fig. 4 shows the evolution of the total magnetic field B to t = 
yj B% + By + B? as a function of time for eight of the stan- 
dard scenarios (-B ga i,o = Bigm,o = 10 -9 G, Fig. 4a & b) and 
for the same scenarios with a lower initial IGM magnetic 
field of Bigm.o = KT 12 G (Fig. 4c & d). We plot the galac- 
tic magnetic field (solid lines) and the IGM magnetic field 
(dotted lines) separately using a threshold of 10 _29 g cm -3 . 
In all of the presented simulations, a slight amplification 
of the total galactic magnetic field by approximately a fac- 
tor of 4 caused by the winding process is visible in the two 
progenitor discs before the first encounter. This behaviour 
was already observed i n simulations of isolated galaxies and 
galactic major mergers (|Kotarba et al.ll2009l . [201ol . l201ll '). All 
of the presented merger simulations show this initial ampli- 
fication with the same order of magnitude. 
During the first encounter, the galactic magnetic field gets 
efficiently amplified within all of the presented simulations, 
whereby the maximum magnetic field strength and the slope 
of the amplification depend on the mass ratio of the progen- 
itor galaxies. Thereby, lower mass ratios lead to a slightly 
higher maximum value of the magnetic field and a steeper 
slope of the amplification, caused by the presumably higher 
turbulence driven by the first encounter. In case of lower 
mass ratios, the maximum value of the galactic magnetic 
field strength is mostly reached after the first encounter. 
However, the maximum value of the magnetic field in sce- 
narios with higher mass ratios is usually not reached until 
the time of the second encounter. 

The spikes in the galactic magnetic field strength after the 
maximum value in the simulations with mass ratios of 2:1 
up to 5:1 (Fig. 4a (inset), 4c) correspond to the second en- 
counters. The second encounter takes place at a later time 
for smaller companion galaxies (cf. section 4.1.). For smaller 
companion galaxies, the second encounter generally results 
in a further increase of the galactic magnetic field strength, 
whereby the maximum value of the magnetic field is reached 
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Figure 3. Evolution of the projection of the mean line-of-sight total magnetic field |B| at eight different time steps (given within the 
panels) for the minor merger scenario M1M4_G6I9. 
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Figure 4. Total magnetic field Btot = y B% + £? 2 + B 2 as a function of time for the eight merger scenarios with the standard 

magnetization of B ga i,o = ^IGM.O = 10" 9 G (upper panels) and for the same eight scenarios but with a smaller initial IGM magnetic 
field of -Bigm.o = lO -12 (lower panels). Galaxies (solid lines) and IGM (dashed lines) are shown separately using a density threshold of 
10" 29 g cm -3 . The evolution of the magnetic field is similar for all merger scenarios, independent of the mass ratios of the progenitor 
discs. However, for smaller companion galaxies, the maximum value of the total magnetic field is slightly lower, and the slope of the 
amplification is slightly flatter compared to the scenarios with larger companion galaxies. 
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Figure 5. Evolution of the total magnetic field Btot = y B| + + _B| as a function of time for different merger scenarios with 

different initial magnetic fields, disc orientations and resolutions, (a) This panel shows the magnetic field evolution for the five merger 
scenarios with the smallest companion galaxies and with the standard magnetization of 10" 9 G. (b) Magnetic field evolution of the 
merger scenarios M1M4 and M1M7, with the highest magnetization of B ga i,o = 10~ 6 G, £?igm,o = 10~ 9 G and also with the smallest 
magnetic field of B ga i q = 10 -9 and £?igm,o = 10" 12 G. (c) Magnetic field evolution as a function of time for the two retrograde scenarios 
M1M3 and M1M7 with the standard initial magnetization of 10 — 9 G together with the corresponding prograde scenarios, (d) Magnetic 
field evolution as a function of time for the scenario M1M3 for the three different resolutions. Galaxies (solid lines) and IGM (dashed 
lines) are again plotted separately using a density threshold of IO -29 g cm~ 3 . 



at the time of the second encounter. However, in most of 
the merger scenarios, the galactic magnetic field decreases 
again after the second encounter. The subsequent encoun- 
ters and the final merger do not lead to any considerable 
further amplification of the galactic magnetic field. At the 
end of the simulations, the galactic magnetic field strengths 
for the presented merger scenarios saturate at a similar value 
of several fiG. The evolution of the IGM magnetic field is 
discussed separately in section 4.2.5. 

Fig. 5 shows the evolution of the total magnetic field Btot = 
\/-B 2 + By + B% as a function of time for different merger 
scenarios with different initial magnetic fields, disc orienta- 
tions and resolutions. Fig. 5a shows this evolution for the five 
merger scenarios with the smallest companion galaxies and 
with the standard initial magnetization of 1CP 9 G. Before 
the first encounter at t = 1.3 Gyr, there is an initial ampli- 
fication of the galactic magnetic field similar to the lower- 
mass-ratio scenarios (cf. Fig. 4, but note the different time 
scales). An exception are the two merger scenarios with the 
smallest companion galaxies M50 and MI 00, which show an 
even stronger initial amplification by approximately a fac- 



tor of 8 due to the winding process. For these simulations 
with very small companion galaxies we used a factor of 10 
higher resolution (see Table 2). The stronger amplification 
can be ascribed to the higher resolution (see chapter 4.2.4 
and 4.4). Within the scenarios involving companion galax- 
ies M10, M15 and M20, the amplification during the first 
encounter results in a galactic magnetic field strength lower 
than 10~ 6 G, which can be ascribed to the lower amount of 
impact energy available for conversion into turbulent energy 
(cf. chapter 4.3.2). The maximum value of the galactic mag- 
netic field strength (1.9 • 10" 6 G, 1.0 ■ 10" 6 G and 0.7 • 10" 6 
G for M1M10, M1M15 and M1M20 scenario, respectively) 
is reached at the time of the second encounter. Similar to 
the scenarios with smaller mass ratios (cf. Fig. 4), subse- 
quent encounters and the final merger do not lead to any 
considerable further amplification of the magnetic field. 

Generally, we find smaller saturation values of the galactic 
magnetic field for mass ratios of 10:1 and higher. Note, how- 
ever, that the merging timescale increases with decreasing 
mass of the companion galaxy, wherefore we were not able to 
simulate the whole merging process for the M1M20 scenario 
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and also the (high resolution) scenarios M1M50 G9I9 and 
M1M100_G9I9. Nevertheless, we can asses a general trend: 
Within the scenarios described above (e.g. the M1M10 and 
M1M15), the galactic magnetic field strength after the fi- 
nal merger does not exceed the maximum value reached af- 
ter the first or second encounter. Therefore we do not ex- 
pect a notable further amplification of the galactic magnetic 
field for the merger scenarios with mass ratios of 20:1 and 
higher. As the high-mass-ratio scenarios M1M50_G9I9 and 
M1M100_G9I9 show only a slight amplification after the 
first encounter (with maximum values within the simulation 
time of8.3-10~ 8 Gand2.6-1(T 8 G for M1M50 and M1M100 
scenario, respectively) , we expect the saturation value of the 
magnetic field to be also of order 10 -8 G. The final magnetic 
field strengths of the merger scenarios with mass ratios of 
10:1 and lower saturate at values of order 10 -6 G, whereas 
scenarios with mass ratios of 15:1 and higher saturate at 
lower values. 

In summary, galaxy mergers with smaller companion galax- 
ies show a slightly lower maximum value of the magnetic 
field strength and a flatter slope of the amplification. At the 
time of the final merger, the saturation values for scenarios 
with a mass ratio up to 10:1 show similar values of order 
/iG, whereas the saturation value for larger mass ratios is 
decreasing with decreasing mass of the companion galaxy. 

4-2.2 Dependence on the initial magnetization 

Fig. 5b shows the total magnetic field as a function of 
time for the merger scenarios M1M4 and M1M7, with the 
highest, i.e. present-day magnetization of -B ga i,o = 10~ 6 
G and -Bigm,o = 10 -9 G. The corresponding simulations 
with the smallest magnetic field of -B ga i,o = 10 -9 G and 
BiGM,o = 10 -12 G are shown for comparison (cf. Fig. 4c & 
d). For the present-day magnetization we find only a small 
peak in the galactic magnetic field evolution at time of the 
first encounter, followed by a slight amplification. The max- 
imum value of the galactic magnetic field strength (5.3 TO -6 
G and 3.8 • 10 -6 G for M1M4 and M1M7 scenario, respec- 
tively) is reached before the second encounter. This maxi- 
mum value is reached at a earlier time and achieves a higher 
value compared to the scenarios with the lower initial mag- 
netic field (4.5 • 10~ 6 G and 2.7- 10 -6 G, respectively). Sub- 
sequently, the magnetic field in the present-day scenarios is 
decreasing. The saturation value within these scenarios is 
comparable to the final value within the simulations with 
lower initial magnetizations. Thus, the initial magnetic field 
strength of the scenarios has indeed an influence on the evo- 
lution of the galactic magnetic field. However, the saturation 
value of the galactic magnetic field is independent from the 
initial galactic magnetic field strength, which is consistent 
with previous studies and will be discussed in more detail in 
section 4.3.2. 



4-2.3 Dependence on the disc orientation 

Fig. 5c shows the magnetic field evolution as a function of 
time for the two retrograde scenarios M1M3 and M1M7 with 
the standard initial mag netization of 10" 9 G (cf. Table 4) 
together with the corresponding prograde scenarios. 
The galactic magnetic field amplification seems to be slightly 



more efficient for retrograde encounters than for prograde 
encounters. This can be ascribed to the presumably more 
efficient turbulence driving during the first encounter, as the 
different rotational velocities within the retrograde scenar- 
ios are expected to serve as additional shearing motions. 
Also, the maximum value of the magnetic field strength is 
higher for the retrograde scenario with the larger compan- 
ion galaxy M3 (7.2 -10 -6 G), compared to the corresponding 
prograde scenario (4.1 ■ 10~ 6 G). The galactic magnetic field 
strength of the retrograde scenario M1M7 is smaller than the 
field strength in the corresponding prograde scenario until 
the second encounter. However, the saturation values of the 
galactic magnetic field within the retrograde scenarios are 
comparable to the values reached in the corresponding pro- 
grade scenarios. Thus, it is interesting to note that prograde 
or retrograde orientation of the galactic disc seems to influ- 
ence the evolution of the galactic magnetic field during the 
collision, but shows only marginal effects on the saturation 
values of the galactic magnetic field. 

4-2-4 Dependence on the resolution 

Fig. 5d shows the magnetic field evolution as a function of 
time for the scenario M1M3 for the three different resolu- 
tions (cf. Table 4). The general amplification behaviour does 
not change significantly with resolution. However, some fea- 
tures in the magnetic field evolution seem to be only resolved 
within the simulation with the highest resolution, particu- 
larly the remarkably steeper increase of the magnetic field 
strength after the first encounter. The simulation with the 
highest resolution shows a higher initial amplification of the 
galactic magnetic field prior to the first encounter. The final 
saturation value of the magnetic field is comparable for all 
resolutions, thus showing that a variation of the resolution 
does not alter the general results. This behaviour is due to 
the better resolution of turbulent structures in simulations 
with higher resolution, while the turbulent energy density 
remains the same for different resolutions. 

4-2.5 IGM magnetic field 

The IGM magnetic field saturates at a value of order nG 
within most of the scenarios. Thus, an initial magnetic field 
of BiGM,o = 10 -9 G does not significantly grow during the 
interaction, whereas an initial field of Bigm,o = 10 -12 G is 
clearly amplified up to the saturation value by three orders 
of magnitude. Within the simulations with an initial galactic 
magnetic field of -B ga i,o = 10~ 6 G, the IGM magnetic field 
amplification at the time of the first encounter is slightly 
more efficient because magnetic energy is additionally trans- 
ported fr om the galaxies into t he IGM by interaction-driven 
outflows (|Kotarba et al.|[201ll ). However, at the end of the 
simulations, we find the same saturation value of the IGM 
magnetic field as within all other scenarios. Moreover, the 
saturation value of the IGM magnetic field is on general in- 
dependent on the masses of the progenitor galaxies up to a 
mass ratio of approximately 10:1. For larger mass ratios, the 
IGM magnetic field saturates at values slightly lower than 
order nG. This indicates that above a ratio of 10:1 the inter- 
action is not violent enough to release a notable amount of 
energy into the IGM and thus to amplify the IGM magnetic 
field. 



12 A. Geng et al. 



Table 5. Impact velocities and energies, maximum magnetic field 
strength reached between first and second encounter 



Collision 


"Gall 


"Gal2 


^impact 


^raax 




[km/s] 


[km / s] 


[10 42 erg] 


[10" 6 G] 


M1M2 


2.92 


3.13 


1.76 


4.6 


M1M3 


2.71 


3.44 


1.49 


4.1 


M1M4 


2.21 


4.12 


1.20 


4.1 


M1M5 


1.95 


4.52 


1.05 


3.4 


M1M6 


1.66 


4.95 


0.91 


3.6 


M1M7 


1.43 


5.28 


0.80 


2.2 


M1M8 


1.29 


5.49 


0.73 


2.1 


M1M9 


1.18 


5.74 


0.68 


1.8 


M1M10 


1.05 


5.93 


0.61 


1.9 


M1M15 


0.84 


6.45 


0.47 


1.0 


M1M20 


0.75 


6.88 


0.39 
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Figure 6. Maximum magnetic field strength B max (reached be- 
tween the first and second encounter within each scenario) as a 
function of the impact energy i?i mp , each normalized to the values 
for the M1M2 merger scenario (cf. Table 5). There is a clear linear 
correlation (green line) between the the maximum magnetic field 
strength and the impact energy of the progenitor galaxies. 
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Figure 7. Evolution of the pressure components Phyd (black 
lines), Ptm-b ( re d lines) and P ma g (green lines) as a function of 
time for the scenarios M1M4 (a) and M1M10 (b) with the stan- 
dard initial magnetization of 10~ 9 G. Pressure values for galaxies 
(solid lines) and IGM (dashed lines) are plotted separately using 
a density threshold of 10 — 29 g cm -3 . At the time of the final 
merger, the magnetic and the turbulent energy densities show 
approximately the same order of magnitude. 



4.3 Self-regulation of the amplification 

As we have shown in section 4.2., encounters of the galaxy 
Ml with a companion galaxy lead to a slightly higher maxi- 
mum value of the total galactic magnetic field the lower the 
mass ratio of the progenitor discs and the larger the compan- 
ion galaxy, respectively. This behaviour can be understood 
in terms of the impact energy of the colliding galaxies. 



4-3.1 Impact energy 

The kinetic energy of the progenitor galaxies, which is re- 
leased during the interaction, is expected to be partly con- 
verted into magnetic energy. Thus, the higher the impact 
energy, the higher the amount of turbulence and the more 
efficient the expected field amplification. The impact energy 
can be estimated on the basis of the masses and the center- 
of-mass-velocities of the progenitor galaxies just before the 
first encounter (this velocity is much less than the initial ve- 
locity of the galaxies since the galactic haloes already over- 
lap, which causes a slowdown of the galaxies). The impact 



energy may be approximated by 

1 2 1 2 

^impact = -m G l«Gl + 2 m G2"G2- (16) 

Center-of-mass-velocities and impact energies for the differ- 
ent merger scenarios are listed in Table 5. As we set the 
galaxies on nearly-parabolic Keplerian two-body orbits, the 
center-of-mass velocity of the companion galaxy increases 
with decreasing mass, whereas the center-of-mass velocity 
of the larger galaxy Ml is slightly decreasing with decreas- 
ing mass of the companion galaxy. As a result, we find the 
impact energy of the M1M5 model to be roughly half the 
energy of the M1M2 model and the energy of the M1M10 
model to be roughly one third of the impact energy of 
our largest merger model M1M2. This is consistent with 
the trend that scenarios with smaller companion galaxies 
show lower maximum values of the galactic magnetic field 
strength (cf. Table 5) and also flatter slopes of the ampli- 
fication during the interaction (cf. section 4.2.). The rea- 
son for this behaviour can most probably be ascribed to 
the lower impact energy which is available for conversion 
into magnetic energy. The resulting correlation between the 



cq 



1.2 

1 

0.8 



°? 0.6 
J 0.4 
0.2 




Magnetic field amplification and X-ray emission in galaxy minor mergers 

(b) 



13 



— 1 1 

galaxies 


- 1 1 1 1 1 

M1M2 G9I9 


IGM 


M1M3~G9I9 




M1M4 G9I9 - 




>u M1M5_G9I9 






■ — ■ 





1.2 
1 

_ 0.8 
« 0.6 
■a 0.4 
0.2 




1 1 

galaxies 


M1M6 G9I9 - 


IGM 


M1M7 G9I9 - 




M1M8 G9I9 - 




M1M9 G9I9 








Figure 8. Mean numerical divergence measure {/i|V • B|/|B|) as a function of time for eight collision simulations with the standard 
initial magnetization of 10~ 9 G (upper panels) and a lower initial IGM magnetic field of Bjqm.o = 10~ 12 G (lower panels). Galaxies 
(solid lines) and IGM (dashed lines) are shown separately using a density threshold of 10 — 29 g cm~ 3 . The numerical divergence measure 
stays below the tolerance value of unity during all of the simulations. 



maximum magnetic field strength and the impact energy is 
shown in Fig. 6. 



4-3.2 Pressures 

The release of kinetic energy during the interaction drives 
turbulence, which in turn results in an amplification of the 
magneti c field by magnetic field line foldin g and shearing 
(see e.g. iBrandenburg fc Subraman ian 2005!, for a review). 
However, if the magnetic energy density reaches the magni- 
tude of the turbulent energy density, the magnetic field am- 
plification caused by the turbulent motion of the gas is sup- 
pressed by the magnetic field itself via the Lorentz force. The 
system therefore tends to maintain a dynamic equilibrium or 
equipartitio n between turbulent and magnetic energy den- 
sity (see e.g. iBe ck 2007; Chv zv et al1l2007l ; lArshakian et all 
120091 ; iKotarba et al.ll2010l . l201ll ). 

In order to study the expected energy equipartition between 
the magnetic and the turbulent energy density (or, equiva- 
lently, magnetic and t urbul ent pressure) in more detail, we 
follow IKotarba et"aH (|201lT ) and choose v IulB (rms velocity 
around the mean velocity inside the smoothing length h) as 



an estimator for the local turbulent velocity 3 . The turbulent 
pressure is then given by P tur b = l/2pv 2 ma . 
Fig. 7 shows the turbulent pressure, the hydrodynamic pres- 
sure -Phyd = l/2pv 2 (with v the total velocity of each parti- 
cle), and the magnetic pressure P mag = B 2 /8iv exemplarily 
for the M1M4 and the M1M10 scenarios with the standard 
magnetization of 10 -9 G. The pressure values for the galax- 
ies (solid lines) and the IGM (dashed lines) are plotted sep- 
arately using a density threshold of 10~ 29 g cm -3 . While the 
hydrodynamic pressure within the galaxies in both scenar- 
ios evolves relatively smoothly and stays in the same range 
of magnitude during the whole interaction (except for the 
peaks indicating the different encounters of the discs), the 
turbulent pressure slightly increases after each encounter. 
During the collision, the galactic turbulent pressure lies be- 
low the galactic hydrodynamic pressure by roughly two or- 
ders of magnitude. The magnetic pressure gets strongly am- 
plified during the interaction (according to the amplification 



»rms is a c omparatively conserv ative estimator for the turbulent 
velocity (see IKotarba et al.ll20TlT ), Hence, it cannot be ruled out 
that the turbulent velocity may be overestimated on small-scales 
and underestimated on larger scales. 
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of the magnetic field itsel f), as has been already show n for 
major galaxy collisions bv lKotarba et all (|2010l . l20Ill ). 
Within both scenarios, the system reaches the equipartition 
level between turbulent and magnetic pressure at about 1 
Gyr after the first encounter. In case of the larger companion 
galaxy in scenario M1M4 (Fig. 7a), the magnetic pressure 
stays slightly above the turbulent pressure until the final 
merger, whereupon an approximate equipartition is reached. 
During the longer period between first and second encounter 
of the scenario M1M10 (Fig. 7b), the magnetic pressure de- 
creases again due to the dilatation of the shocked region, 
which causes also a dilatation of the magnetic field energy. 
However, in the subsequent evolution the magnetic pressure 
gets amplified again with each further encounter, whereas 
each amplification period is followed by a period of decreas- 
ing magnetic field until the next encounter takes place. At 
the time of the final merger, the magnetic and the turbu- 
lent energy densities show approximately the same order of 
magnitude. As the equipartition value within each model de- 
pends on the energy supplied to the system (impact energy), 
we generally find higher equipartition values for lower mass 
ratios of the galaxies. 

Within each of the scenarios shown in Fig. 7, the IGM shows 
a slight increase in the hydrodynamic pressure at the begin- 
ning of the simulation, followed by a relatively smooth evolu- 
tion. The IGM turbulent pressure clearly shows some peaks, 
corresponding to the encounters. The IGM turbulent pres- 
sure lies below the IGM hydrodynamic pressure by roughly 
two orders of magnitude during the whole simulation, which 
is a comparable separation as for the corresponding galactic 
pressure components. The IGM magnetic pressure within 
each scenario gets amplified during the first encounter. In 
the subsequent evolution, a loose equipartition between the 
IGM turbulent and magnetic pressure is maintained until 
the end of the simulations. 

In summary, within all of our merger simulations, we see a 
general trend: magnetic and turbulent pressures of galaxies 
and IGM tend to reach equipartition, which is consistent 
with studies of major galaxy collisions (|Kotarba et al.|[2010l . 
l201ll) and in good agr eement with theoretical expectations 
l|Arshakian et alj|200a ). 

4.4 Numerical divergence and resolution study 

In SPMHD simulations the numerical divergence measure 
(h\ V-B]/|B]) is generally regarded as a measure of reliability 
of the simulations. It is important to note that the numerical 
divergence is not a physical dive rgence caused by magnetic 
monopoles. lKotarba et al] (|2009l ) have demonstated that the 
numerical divergence in simulations of isolated galaxies per- 
formed with Euler Potentials (which are free of a physical 
divergence by construction) can still reach values up to order 
of unity. Hence, the value of the numerical divergence mea- 
sure in SPMHD simulations with direct magnetic field im- 
plementation should not exceed unity in order to guarantee 
the rel iability of the s i mulat ion r esults. As already po inted 
out by IKotarba et all (|201fj| ) and iBiirzle et all (|2011al ). the 
SPH estimator for the divergence (V • B)i (with i being the 
index of the particle) accounts for weighted magnetic field 
differences of neighbouring particles within the smoothing 
length hi, and is therefore primarily dependent on the irreg- 
ularity of the magnetic field within the kernel. In order to 
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Figure 9. Mean numerical divergence measure (h\ V • B|/|B|) as 
a function of time for the M1M3 scenario with the standard mag- 
netization of 10 — 9 G, using three different resolutions. Galaxies 
(solid lines) and IGM (dashed lines) are shown separately using 
a density threshold of 10~ 29 g cm~ 3 . The divergence measure is 
higher for lower resolutions. 



avoid instabilities due to the numerical divergence errors, a 
subtraction of any non-vanishing divergence from the mo- 
mentum equation (see chapter 2.2) is adopted for the simu- 
lations presented in this paper. Also, the adopted IGM helps 
to increase the numerical stability, as the magnetic field is 
not dropping to zero at the edges of the galaxies, which could 
lead to incorrect calculations of the numerical divergence in 
these regions. 

Fig. 8 shows the mean numerical divergence measure (ft|V ■ 
B]/|B|) as a function of time for eight different merger sce- 
narios with the standard initial magnetization of 10 -9 G 
(Fig. 8a & b) and a lower initial IGM magnetic field of 
Bigm,o = 10~ 12 G (Fig. 8c & d). The values of the numerical 
divergence stay below unity within all of the simulations and 
approach unity only during the first encounter. On general, 
the divergence measure is higher for smaller initial magneti- 
zations, in agreement with previous studies. This behaviour 
follows from the fact that in case of weaker magnetic fields 
the encounters exert stronger effects on the irregular mo- 
tions of the particles, wh ich in turn are respo nsible for the 
numerical divergence (cf. IKotarba et ai1l201fj| ). 
Fig. 9 shows the mean numerical divergence measure {h\X7 ■ 
B|/|B|) as a function of time for the three different resolu- 
tions of scenario M1M3 with the standard magnetization of 
I0~ 9 G. For the high resolution (green lines), the numeri- 
cal divergence measure is by a factor of two lower than for 
the medium resolution (blue lines), and for the low reso- 
lution (red lines), it is by a factor of two higher than for 
the medium resolution. This can be explained considering 
the smoothing length hi, which is smaller for a higher reso- 
lution. For a smaller smoothing length, we expect less sub- 
resolution tangling of the magnetic field, resulting in a lower 
numerical divergence measure for a higher resolution. 
On the other hand, the saturation value of the magnetic 
field does on general not depend on the resolution (Fig. 5d), 
showing that even a by a factor of two higher numerical 
divergence would not alter the general results. Hence, we 
conclude that the presented medium-resolution simulations 
are numerically reliable. 
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Figure 10. Evolution of the projection of the mean line-of-sight temperature at the same time steps and for the same scenario as in 
Fig. 1. 



4.5 Temperature evolution and X-ray emission 



The X-ray emission during a galactic encounter is expected 
to depend on the orbit, the disc orientatio n, the masses an d 
the gas fraction of the progenitor discs (e.g. lCox et al.ll2006j ). 
The main source of X-ray emission is hot gas, which is heated 
by shocks accompanying the interaction. As the gas density 
within the IGM and the strength of the magnetic field have 
an impact on the gas behaviour during a galactic interaction, 
they may be expected to also alter the X-ray emission by a 
non-negligible order of magnitude. 



Table 6. Maximum and final values of IGM and total tempera- 
ture and total X-ray luminosity corresponding to Fig. Ill I 



As shown in iKotarba et all (|201lf ), the propagation veloc- 
ity and thus the Mach numbers of shocks driven by galac- 
tic collisions are higher for stronger initial magnetic fields. 
This is due to the magnetic pressure, which is supporting 
the shocks and galactic outflows. Also, the higher this ad- 
ditional pressure component, the less the gas can compress, 
thus leading to lower gas densities. Higher Mach numbers 
are in turn related to the temperature of the IGM behind 
the shock fronts (Rankine-Hugoniot shock jump conditions). 
Additionally, the region of gas which is heated by the shock is 
larger for a higher shock propagation velocity and a stronger 
galactic outflow, respectively. As a consequence, the X-ray 
emission may be expec ted to increase with increasing mag- 
netic field strength (cf. IKotarba et afl|201ll ). 
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Tmax 
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[10 5 K] 


[10 s K] 


[erg s _1 ] 


[erg s _1 ] 
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M1M4_ 


G6I9 
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38.57 


36.56 


M1M4_ 


G9I9 


6.84 


4.29 


38.52 


36.78 


M1M4_ 


G0I0 


5.88 


4.56 


38.49 
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M1M7_ 


G6I9 
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38.55 
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M1M7_ 


G9I9 


6.12 
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38.50 


36.42 


M1M7_ 


G0I0 
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TOTAL 
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4.53 
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G0I0 
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5.53 
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M1M7_ 


G6I9 


6.99 


4.12 


38.78 


36.68 


M1M7_ 


G9I9 


6.49 


4.69 


38.74 


36.73 


M1M7 


G0I0 


6.36 


5.40 


38.80 


38.15 
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We want to investigate the dependence of the X-ray emis- 
sion on the magnetic field strength in the presence of 
an ambient IGM (thus expanding previous stud ies, e.g. 
ICox et all |2006| ; ISinha fc Hollev-Bockelmannl |2009h . There- 
fore, we have calculat ed the bolometric X-r ay luminosity fol- 
lowing the method of iNavarro" et all (| 19951 ). This method is 
based on the assumption that the main X-ray mechanism is 
thermal bremsstrahlung, in agreement with the applied zero- 
metallicity cooling function (cf. section 2.3). The bolometric 
X-ray luminosity is calculated according to jNavarro et al.l 



= 1.2 ■ 10" 



E 



k B Tj 
keV 



1/2 



(17) 

with mass m g a S ,i, density pi and temperature T; of the i- 
th gas particle in cgs-units, respectively Only fully ionized 
particles should be considered when calculating the lumi- 
nosity Therefore, we exclude contributions of particles with 
temperatures lo wer than 10 5 2 K and densities higher than 
O.OIMqpc" 3 (cf. lCox et afll2006f ). 
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Figure 11. Mean temperature (left panels a&c) and mean bolometric X-ray luminosity (calculated according to Eq. (17), right panels 
b&d) of the fully ionized particles as a function of time for the merger scenarios M1M4 {solid lines) and M1M7 {dotted lines), and with 
different initial magnetizations of -B ga i,o = 10 — 6 G, £?igm,o = 10 — 9 G (black lines), S ga i,o = B\GM,0 = 10 _9 G (red lines) and excluding 
magnetic fields (green lines). The upper panels show the temperature and luminosity of the IGM (applying the same threshold as in Fig. 
4-7), and the lower panels show the total (galaxies plus IGM) quantities. The higher the initial magnetization, the more efficient the 
increase of the temperature and thus the X-ray luminosity within the IGM after the first encounter at t m 1.3 Gyr. However, the presence 
of magnetic fields does not enhance the total X-ray luminosity. After the first encounter, the total X-ray luminosities decrease faster and 
result in a lower final value of the luminosity for higher initial magnetizations. This may be explained by the stronger interaction-driven 
outflows within the simulations with higher initial magnetizations, which lead to a faster dilution of the IGM gas and thus a lower gas 
density within the X-ray emitting regions. 



4-5.1 Temperature evolution 

Fig. [10] shows the morphological evolution of the projection 
of the mean line-of-sight temperature for the present-day 
scenario M1M4_G6I9, whereby the eight panels display a 
sequence of increasing time. Before the first encounter at 
about 1.3 Gyr, the galaxies are surrounded by a small region 
of hot gas due to the initial accretion of IGM gas onto the 
galaxies, whereas the gas within the galaxies is cooler (Fig. 
I lOt left upper panel). During the first encounter, shocks and 
outflows are driven into the IGM, thus heating the IGM ga s 
(Fig. 1101 second upper panel) (cf. also lKotarba et al]|201ll ). 
After the first encounter, the temperature of the gas sur- 
rounding the galaxies decreases again due to the dilatation 
of the shocked IGM (Fig. 1101 left lower panel). Subsequent 
encounters result in a further shock-heating of the IGM gas. 



After the final merger (f» 4 Gyr), the IGM gas slowly cools 
down again. 

Fig. [TJJ shows the temperature (left panels) and the bolo- 
metric X-ray luminosity (right panels) of the fully ionized 
particles, respectively, as a function of time for the merger 
scenarios M1M4 (solid lines) and M1M7 (dotted lines), and 
with different initial magnetizations of -B ga i,o = 10~ 6 G, 
Bigm.o = 10~ 9 G (black lines), B ga i,o = Bigm.o = 10~ 9 G 
(red lines) and excluding magnetic fields (green lines). The 
upper panels show the temperature and luminosity of the 
IGM, and the lower panels show the total (galaxies plus 
IGM) quantities. Table 6 summarizes the maximum tem- 
peratures (T max ) and X-ray luminosities (Lx max ) reached 
during the interaction, as well as the final values (Tfi na i and 
Lxfi na i) reached at t — 16 Gyr (cf. Fig. [TJJ. 

Fig. 1 1 lb , shows the mean IGM temperature as a function of 
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Figure 12. Evolution of the projection of the mean line-of-sight X-ray luminosity at the same time steps and for the same scenario as 
in Fig. 1. 



time, applying the same density threshold as before in Figs. 
4-7. We have set the initial IGM temperature to the virial 
temperature of the larger galaxy Ml (cf. section 3.4), thus 
overestimating the equilibrium temperature of the merging 
system. Thus, the IGM gas tends to cool in the very begin- 
ning of the simulations. At the time of the first encounter 
(t « 1.3 Gyr), interaction-driven shocks reheat the IGM gas, 
resulting in temperatures considerably higher than the ini- 
tial temperature. Thereby, the increase of the IGM temper- 
ature is more efficient for stronger initial magnetizations, 
in ag reement with the previous studies (cf. iKotarba et al.l 
l201lf ). reaching higher maximum temperatures of the IGM 
for higher initial magnetic field strengths (cf. Table 6). Fur- 
thermore, the maximum temperature is higher for larger 
companion galaxies, i.e. for higher impact energies (cf. Ta- 
ble 5). After the final merger (t « 4 Gyr for the M1M4 and 
t ~ 6.5 Gyr for the M1M7 scenario) the IGM temperature 
decreases. The final IGM temperatures lie between 4 ■ 10 5 K 
and 4.7 • 10 5 K within all of the scenarios shown in Fig. 11 
(see also Table 6). 

Fig. II lb shows the mean total temperature as a function 
of time. The behaviour of the total temperature is gener- 
ally very similar to the behaviour of the IGM temperature, 
because the interaction-driven shocks propagate mainly into 
the IGM and deposit their energy there. The maximum total 
temperatures reached during the first encounter are slightly 
higher than the maximum IGM temperatures. However, the 
spread in the final total temperatures is considerably larger 
compared to the IGM values (see Table 6). 



4-5.2 Evolution of bolometric X-ray emission 

Fig. [T2] shows the morphological evolution of the projection 
of the mean line-of-sight X-ray emission for the present-day 
scenario M1M4_G6I9, whereby the eight panels display the 
same sequence of increasing time as in Fig. 1. Before the 
first encounter (upper left panel), the X-ray luminosity is 
relatively high within the galaxies (high densities) and the 
surrounding region (higher temperatures due to the initial 
accretion of the IGM gas, cf. the previous section). Until 



the final merger (t » 4 Gyr for the M1M4 and t « 6.5 
Gyr for the M1M7 scenario), the X-ray luminosity evolves 
according to the temperature evolution, i.e. decreases until 
the first encounter due to the initial cooling of the IGM gas, 
increases during the first encounter due to the shock-heating 
of the IGM, and subsequently decreases due to gas cooling 
and dilatation of the shock-heated regions. 

Fig. II 1 b shows the bolometric luminosity as a function of 
time for the IGM. The correlation between luminosity and 
temperature (Fig. 11a) until the final merger is clearly visi- 
ble. Particularly, the luminosity increases considerably dur- 
ing the first and second encounters (tfi rst ~ 1.3 Gyr and 
tsecond ~ 3.3 Gyr for the M1M4 and isecond ~ 4.3 Gyr for 
the M1M7 scenario, respectively, inlay in Fig. 10b), whereby 
higher initial magnetic field strengths result in a faster and 
stronger increase of the IGM X-ray luminosity directly after 
the first encounter and also higher maximum values. How- 
ever, after the final merger, the IGM X-ray luminosity de- 
creases even more than the IGM temperature. This is due to 
the linear dependence of the luminosity on the density (cf. 
Eq. 17) and the dilution of the shock-heated IGM gas (not 
shown). The final values of the IGM X-ray luminosity lie in 
the range between 36.0 erg s _1 and 37.5 erg s _1 for all of 
the merger scenarios. Most interestingly, the final luminosi- 
ties are higher for lower initial magnetizations. Additionally, 
they are also higher for larger companion galaxies, i.e. higher 
impact energies (see Table 5). We will discuss this behaviour 
in more detail in the next section. 

Fig. Illt i shows the mean total bolometric luminosity of both 
the IGM and the galaxies as a function of time. As within 
the IGM, the behaviour of the total X-ray luminosity follows 
the behaviour of the total temperature until the final merger 
(thus also reaching higher maximum X-ray values compared 
to the IGM maximum values) and subsequently decreases 
faster than the temperature due to the additional dilatation 
of the IGM gas. The final values of the total X-ray luminosity 
(t ps 16 Gyr) lie in the range between 36.7 erg s _ and 38.3 
erg s~ , i.e. are slightly higher as for the IGM only. As before, 
the final luminosities are higher for lower magnetizations and 
larger companion galaxies, i.e. higher impact energies. 
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4-5.3 Discussion of X-ray emission 

The presence of magnetic fields shows a clear influence on 
the evolution of the total X-ray luminosity in the presented 
simulations (Fig. 11). On the one hand, interaction-driven 
shocks have been shown to gain hig her Mach numbers i n 
the presence of a strong magnetic field (|Kotarba et al.|[201ll ). 
thus leading to a more efficient shock-heating of the IGM gas 
and consequently a higher X-ray luminosity within the IGM. 
In this sense, this is a positive feedback of the magnetic field 
on the X-ray luminosity. However, the X-ray emission of the 
galaxies themselves outweights the emission of the IGM at 
the times of the different encounters (see the inlays in Fig. 
11 b & d, particularly the clear peaks in Fig. lid). On the 
other hand, there is also a negative magnetic feedback: The 
stronger the magnetization of the system, the higher the 
magnetic pressure pushing the gas apart, and thus the more 
efficient the dilution of the shock-heated gas. Given that 
the X-ray emission scales with the root of the temperature 
but squared with the gas density (Eq. 17), this negative 
feedback soon overtakes the positive, finally resulting in a 
clearly lower final X-ray luminosity of the system in the 
presence of stronger magnetic fields (Fig. 11 and Table 6). 
This is contrary to what was expected by iKotarba et al.l 
(|201ll ). 

In general, we find relatively low final values of the total 
X-ray luminosity (38.15 — 38.33 erg s _1 in our simulations 
without magnetic field s) compared to the final X-ra y lumi - 
nosity values found in ISinha fe Hollev-Bockelmannl (|2009f ) 
(« 38.5 — 41.5 erg s _1 , depending on the gas content and 
the mass of the progenitor galaxies, cf. their Fig. 6) and 
ICox et all (|200fJ ) (« 39-42.5 erg s~\ depending on the pro- 
genitor mass, cf. their Fig. 11). This is most probably due to 
the lack of an accreting black hole in our models, and that 
zero-metallicity cooling is used for the X-ray co mputation 
in con trast to the metal- line cooling 4 applied bv lCox et alj 
(|2006l ). Furthermore, as the overall X-ray emission depends 
also on the gas con tent, the orbit and the disc orientation 
(cf. ICox et alj|2006h . the initial setup of the scenarios has 
an additional influence on th e X-ray luminosity. 
ISinha fc Hollev-Bockelmannl (|2009T l simulated interacting 
galaxies including hot gas within the galactic haloes, but 
no gaseous disc. This is a clear difference in the setup com- 
pared to our gaseous discs. However, despite this difference, 
we find the X-ray luminosity to increase and reach its maxi- 
mum prior to the final merge r . simi lar to what was found by 
ISinha fc Hollev-Bockelmannl (|2009r ). Furthermore, this evo- 
lution of the X-ray emission corresponds also to observa- 
tional results, which show the X-ray luminosity to increase 
during t he merger and peak abo ut 300 Myr before the final 
merger dBrassington et al]|2007l ). 

ICox et ail (|2006t ) considered the X-ray emission generated 
during interactions of galaxies including gaseous discs, but 



4 Because our intention is to investigate the influence of mag- 
netic fields on L x boloi we are focusin g on the more co nservative 
zero-metallicity method. As shown bv lCox et al.l J2006f) . the pro- 
gression of the X-ray emission is very similar for both methods, 
but the metallicity-dependent cooling leads to a higher X-ray lu- 
minosity during the collisions (see their Fig. 2). Thus, the zero- 
metallicity X-ray luminosity may be regarded as a lower limit 
iCox et alj|2006h . 



without surrounding IGM gas. In good agreement with their 
results, we generally find a higher X-ray luminosity for larger 
progenitor galaxies (cf. our Fig. 11 and their Fig. 11), due 
to the aspired hydrostatic equilibrium after the final merger. 
On the other hand, we find the maximum X-ray luminosity 
after the first or the second encounter , not a t the time of 
the final merger as found bv lCox et all (|2006t ). 
Furthermore, our minor merger simulations show much 
longer time periods between the encounters, thus enabling 
cooling of the gas. This is different to the major merger stud- 
ies with re l ativel y short but violent cooling periods found by 
ICox et ail (|2006T ). Altogether, we find at least o ne order of 
magni tude lower X-ray luminosities compared to lCox et al.l 
(|2006l ) and a faster decrease in luminosity after the final 
merger. Again, this difference may result from the lack of 
black holes and metal-line cooling in our simul a tions , and 
the lack of IGM gas in the studies bv lCox et all (|2006t ). 
In summary, we may underestimate the total X-ray luminos- 
ity within our simulations, but our results generally agree 
well with previous studies and observations. Therefore, we 
are confident that our conclusions on the impact of the mag- 
netic field on the ability of a system to emit in X-ray are 
reliable. 



5 CONCLUSION AND DISCUSSION 

We have presented a series of 32 galaxy minor mergers in- 
cluding galactic magnetic fields and a magnetized ambient 
IGM. We have investigated the evolution of the galactic and 
IGM magnetic field with respect to different mass ratios, ini- 
tial magnetic field strengths, disc orientations and numerical 
resolutions. The main results can be summarized as follows: 

• The smaller the companion galaxy, i.e. higher the mass 
ratio, the lower the maximum galactic magnetic field and 
the flatter the slope of the amplification. The magnetic field 
for mass ratios up to 10:1 saturates at values of order /xG, 
whereas scenarios with smaller companion galaxies saturate 
at values lower than fiG, suggesting that the impact energies 
for mass ratios above 10:1 are not high enough to provide 
enough energy for conversion into magnetic energy. 

• The saturation value of the galactic magnetic fields is 
independent on the initial galactic magnetic field strength, 
in good ag reement with previous studies (|Kotarba et al.l 
l2010l . l201ll ). 

• The disc orientation (prograde or retrograde) slightly 
changes the evolution of the galactic magnetic field during 
the encounter, but has only a minor effect on the saturation 
value of the galactic magnetic field. 

• Some features of the magnetic field evolution seem to be 
only resolved in simulation with the highest numerical res- 
olution. However, the saturation value of the magnetic field 
is comparable for the three investigated resolutions. Most 
probably, this is due to the better treatment of small-scale 
turbulence in high-resolution simulations, which, however, 
does not significantly increase the total turbulent energy of 
the system. 

• The IGM magnetic field saturates at a value of order 
10 _9 G within most of the scenarios. For mass ratios higher 
than 10:1, the IGM saturation value is slightly lower, which 
can again be traced back to the lower impact energy of these 
interactions. 
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Summing up, minor mergers up to a mass ratio of 10:1 are 
able to amplify a small (10~ 9 G) galactic magnetic field 
up to fiG order, a val ue which is in good agreement with 
observations (see e.g. iBeck et all Il996l ; IChvzv et alj 120071 ; 
IVollmer et aT]|2010l ). Also, the saturation value of the IGM 
magnetic field of ~ 10 -9 G is in very good agreement with 
the upper li mit estimates of 10~ 9 — 10 -8 G derived from 
observations (|Kronberg et al.ll2008l ). 

Furthermore, the maximum values of the galactic magnetic 
field reached during the interactions are higher and the slope 
of the amplification is steeper for lower mass ratios. This is 
reasonable because the main source for the magnetic field 
amplification is the impact energy released during the inter- 
action. This impact energy is in turn higher for lower mass 
ratios (section 4.3.1 and Table 5). 

Moreover, we find that the magnetic pressures associated 
with the IGM and galactic magnetic fields saturate at the 
equipartition level between turbulent and magnetic pressure 
(section 4.3.2), w hich is consistent w i th st udies of major 
galaxy collisions (|Kotarba et al.l |2010| . l201lh and theoreti- 
cal considerations. As the energy equipartition within each 
scenario depends on the impact energy supplied to the sys- 
tem, the equipartition values are on general higher for larger 
companion galaxies, i.e. lower mass ratios. 
The values of the numerical divergence measure /i|V-B|/|B| 
(which are generally non-zero within SPMHD simulations) 
have been shown to lie below the tolerance-threshold of unity 
within all of our simulations (section 4.4). Thus, our sim - 
ulations are numerically reliable (cf. iKotarba et al.ll2010l ). 
Moreover, our resolution study showed that the numerical 
divergence measure is lower for a higher numerical resolu- 
tion, but simultaneously the saturation value of the mag- 
netic field is generally independent of the resolution (section 
4.2.4). Thus, our conclusions are robust with respect to the 
numerical divergence and resolution. 

Finally, we presented a detailed analysis of the evolution of 
the temperature and the X-ray luminosity during the sim- 
ulated minor merger scenarios (section 4.5). We find that a 
magnetic field has both, a positive and a negative feedback 
on the X-ray luminosity. On the one hand, a higher mag- 
netic field results in higher Mach-numbers of the interaction- 
driven sh ocks and thus a mor e efficient shock-heating of the 
IGM (cf. IKotarba et alJlamh . Consequently, the X-ray lu- 
minosity is higher within recently shock-heated regions. On 
the other hand, a higher magnetic pressure leads also to a 
more efficient dilution of the shock-heated regions, result- 
ing in lower gas densites and hence lower X-ray luminosi- 
ties. Given the different dependence of the X-ray luminos- 
ity on the temperature and the gas density (Eq. 17), the 
negative feedbeck is more important in the long run. Con- 
sequently, we find lower final X-ray luminosities for higher 
initial magnetization s, in contrast to what was expected by 
IKotarba et"afl (|201lh . 

As galaxy mergers are believed to be an essential part of hi- 
erarchical growth of structure in the universe, these events 
may be expected to provide a non-negligible contribution 
to the amplification of magnetic fields on galactic scales 
and, given the simultaneous magnetization of the ambient 
IGM, also beyond. Commonly, the presence of galactic mag- 
netic fields is explained by the action of the galactic dy- 
namo (e.g. lFerriere|[l992l ;l Gressel et al.l [20081 ; iHanasz et al.l 
l2009h . However, this hypothesis is challenged by recent ob- 



servations of strong (fiG) magnetic fields at high redshifts 
llBernet et al. 1 120081 1, because the galactic dynamo can am- 
plify magnetic fields only on a timescale of several Gyr. Also, 
it has been shown that that dynamos cannot work efficiently 
within d warf galaxies due to the lack o f global differential 
rotati on (|Gressel et al.l [20081 ) . However. iDubois fc Tevssierl 
presented a possible machanism for the magnetiza- 
tion of the early universe. Within their "Cosmic Dynamo" 
the universe is magnetized by the combined action of grav- 
itational instabilities forming dwarf galaxies, a galactic dy- 
namo working within these newly formed dwarf galaxies, 
and galactic winds expelled into the IGM at times of violent 
SF activity. 

Nevertheless, galactic interactions provide a promising al- 
ternative or complement to the scenarios proposed so far: 
As galaxy mergers, and especially minor mergers, are be- 
lieved to have been much more frequent in earlier times of 
the universe, it is likely that the presented interaction-driven 
amplification also provided for a significant contribution to 
the amplification of the galactic and intergalactic magnetic 
fields on short timescales. Hence, for future studies it would 
be interesting to focus on mergers of dwarf galaxies to gain 
a more complete picture of the evolution of magnetic fields 
within mergers. Furthermore, studies of the magnetic field 
evolution in the early universe in a cosmological context and 
within the formation of galaxies would be very worthwhile 
for a better understanding of the overall process of the mag- 
netic field amplification caused by galaxy formation and in- 
teraction events in the history of the universe. 
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